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Abstract : Two similar Minkowskian diffusions have been considered, on one hand 
by Barbachoux, Dcbbasch, Mahk and Rivet ([BDRl], [BDR2], [BDR3], [DMR], [DR]), 
and on the other liand by Dunkel and Hiinggi ([DHl], [DH2]). We address here two 
questions, asked in [DR] and in ([DHl], [DH2]) respectively, about the asymptotic 
behaviour of such diffusions. More generally, we establish a central limit theorem 
for a class of Minkowskian diffusions, to which the two above ones belong. As a 
consequence, we correct a partially wrong guess in [DHl]. 

1 Introduction 

Debbasch, Malik and Rivet introduced in [DMR] a relativistic diffusion in Minkowski space, 
they called Relativistic Ornstein-Uhlenbeck Process (ROUP), to describe the motion of a point 
particle surrounded by a heat bath, or relativistic fluid, with respect to the rest-frame of the fluid, 
in which the particle diffuses. This ROUP was then studied in [BDRl], [BDR2], [BDR3], [DR], 
and extended to the curved case in [D]. Then Dunkel and Hiinggi introduced and discussed in 
[DHl], [DH2] a similar process, in Minkowski space, they called relativistic Brownian Motion. 

Note that independently, a Relativistic Diffusion on any Lorentz manifold was defined in [FLJ], 
as the only diffusion whose law possesses the relativistic invariance under the whole isometry group 
of the manifold. Accordingly, the particular case of the Schwarzschild-Kruskal-Szekeres manifold 
was studied. The case of Godel's universe was recently studied in [F]. 

In [DR], Debbasch and Rivet argue qualitatively that the so-called "hydro dynamical limit" of 
their ROUP should behave in a Brownian way. They stress that a mathematical rigourous proof 
remains needed, to confirm such not much intuitive statement. 

In [DHl] and [DH2], Dunkel and Hanggi ask the question of the asymptotic behaviour of 
the variance, or "mean square displacement", of their diffusion. Indeed, comparing to the non- 
relativistic case, and after numerical computations, they guess that this variance, normalised by 
time, should converge, to some constant for which they conjecture an empirical formula. 

We answer here these two questions, asked by Debbasch and Rivet in [DR], and by Dunkel 
and Hanggi in [DHl], [DH2], and indeed a more general one. We establish in fact rigourously 
in this article, a central limit theorem for a class of Minkowskian diffusions, to which the two 
above mentioned ones, ROUP and Dunkel-Hanggi (DH) diffusion, belong. As a consequence of 
our main result, we establish, for this whole class, the convergence of the normalised variance, 
guessed (for their particular case) in [DHl] and [DH2]. Getting the exact expression for this 
limiting variance, and particularising to the (DH) diffusion, we can then invalidate and correct 
the wrong conjecture made in [DHl] on its expression and asymptotic behaviour (as the noise 
parameter goes to infinity). 

To summarise the content, we begin by describing in Section 2 below the class of Minkowskian 
diffusions we consider, which contains both ROUP and DH diffusions as particular cases. 

Then in Section 3, we present our study, leading to the following main result : 
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Theorem 1 Let (x(,pt)(>o = {xl,pl)i<i<d,t>o 6e a M'^ x W^-valued diffusion solving the sto- 
chastic differential system : for 1 < i < d , 

' dxi = /(|pt|)pj(it 

^""^ I dp| = -6(|pi|)pjdt + a(|pt|)(/3[l + r?(|pj|)2])"'^'[dT^/ + 7?(|p,|)(pj/|p,0^ ' 

Then, under some hypotheses (TC) on the continuous functions f, b, a , rj , the law of the process 
Xat)^^g converges, as t oo , to the law of {T,pBa)^yQ, in C(M"'",M'^) endowed with 
the topology of uniform convergence on compact sets of M+ . Here W and B are standard d- 
dimensional Brownian motions, w is a standard real Brownian motion, independent ofW, f3 > 
is an inverse noise or heat parameter, and is a constant, displayed in Proposition 2 below. 

The fohowing is then deduced. Recah that the symbol E stands for expected or mean value, 
with respect to the underlying probability measure (or distribution) ¥ (governing the given 
Brownian motions W,w). 

Corollary 2 Under the same hypotheses as in the above theorem, from any starting point, the 
normalised variance (mean square displacement) t~^ E [jx^p] goes, as t — > oo , towards d x . 

In Section 4, we study the behaviours of the limiting variance S| , as the inverse noise para- 
meter 13 goes to or to oo , and we also support our result by numerical simulations. Focussing 
on the particular case of the DH diffusion, for d = 1 , we get the following, which, though confir- 
ming a non-classical variance behaviour, shows up a behaviour near which differs from the one 
implied by the wrong guess made in [DHl] about the expression of S| . 

Proposition 4 Consider the DH case, for d = 1 , as in [DHl]. Then, we have 

So ~ 2//3 as [3 y oo ; and, as (] \0 : — — for some explicit constant A > . 

^ ^ log(l//3) 

Finally, we detail in Sections 5 and 6 two somewhat involved proofs. We thank Reinhard 
Schafke for his kind and decisive help for the proof of Sections 5.1 and 5.2. 



2 A class of Minkowskian diffusions 

Let K^'"^, where d>l is an integer, denote the usual Minkowski space of special relativity. 
In its canonical basis, denote by x = (x^) = {x^ , x*) = (x*^, x) the coordinates of the generic 

point, with greek indices running 0, and latin indices running l,..,(i. The Minkowskian 

d 

pseudo-metric is given by : ds"^ = jdx'^p — ^ |(ix*|^. 

i=l 

The world line of a particle having mass m is a timelike path in M^'*^, which we can always 
parametrize by its arc-length, or proper time s . So the moves of such particle are described by 
a path s I— > (xs), having momentum p = (ps) given by : 

P = (P^) = (P°, P') = (P°, P) , where := m ^ , 

0,2 ^,..,2^^2_ 



and satisfying : Ip'^P — \p' 



2 



We shall consider here world lines of type {t, x(t))t>o i take m = 1. Introducing the 
velocity v = {v^, . . . , v'^) and polar coordinates (r, Q) by setting : 

v'--=^, r:=\p\={Y,\pr] and @ := ^ =: {e\ . . . , O'^) e 
we get at once : 

^0 



,i=l 



p" = ^ = Vl + r2= (l-lvp) and p = {1, . 

as 

Thus, a full space-time trajectory 

p{t))t>o = {t, x(t), p°{t), p{t))t>0 

is determined by the mere knowledge of its spacial component (x(t), p(t)). 

We can therefore, from now on, focus on spacial trajectories t (xt, pt) G 1^*^ x 

The Minkowskian diffusions we consider here are associated as above, to Euclidian diffusions 
t I— > (xf, Pt) = (2;J,pJ)i<i<(i , which are the solution to a stochastic differential system of the 
following type : 

dxl = f{rt)p\dt 

/ N -1/2 , for 1 < z < d, 

dp\ = - b{rt)p\ dt + a{rt) [15 [1 + r/(rj)2]j [dWt + 7?(rt) 61 dwt] 

where W := (W^^.-.^W^) denotes a standard d-dimensional Euclidian Brownian motion, w 
denotes a standard real Brownian motion, independent of W, /? > is an inverse noise or heat 
parameter, and the real functions /, 5, <t , rj are continuous on ]R+ and satisfy the following 
hypotheses, for some fixed e > : 

iji) cr > e on M+ ; gir) := ^ ^^y} > e for large r ; lim ^ fif) = for some e' < 4^ . 

(7 (r) r—>oo 



Of course, in the particular case of constant functions /, b, a, and r] = 0, the process (x^) 
is an integrated Ornstein-Uhlenbeck process. The process considered by Debbasch, Malik and 
Rivet ([BDRl], [BDR2], [BDR3], [DMR], [DR]), they call Relativistic Ornstein-Uhlenbeck Process 
(ROUP), corresponds to : 

{ROUP) /(r) = &(r) = (l + r2)-V2^ ^(^)^^^ ^^0, g{r) = r {1 + r^)'^/^ , 
and the relativistic process considered by Dunkel and Hanggi ([DHl], [DH2]) corresponds to : 

(DH) f{r) = {1 + r^)-^/^ , b{r) = 1, a{r) = \/2y^Y^ , 7]{r) = r , g{r) = r {1 + r^)'^/^ . 

These processes are intended to describe the motion of a point particle surrounded by a heat 
bath, or relativistic fluid, with respect to the rest-frame of the fluid, in which the particle diffuses. 
The Minkowskian diffusion (xf, pf) solving the stochastic differential system (★) is isotropic pre- 
cisely when rj = (for d > 2; when d = 1, r] does not matter). If rj ^ , the momentum (pt) 
undergoes a radial drift. 
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In ([DR], Section 4), Debbasch and Rivet argue heuristically that the so-cahed "hydro dynami- 
cal Hmit" of their ROUP should behave in a Brownian way, and ask the question of a mathematical 
proof confirming such not much intuitive statement. 

In ([DHl], [DH2]), Dunkel and Hanggi ask the question of the convergence, as t goes to 
infinity, of the normalised variance (or mean square displacement) : 



:= r^E ( ^ 



E 



■,1=1 



We shall answer these two questions, by means of the more general one we address, which is 
the asymptotic behaviour, as t ^ oo , of the process : 



>0 



at } = t 

a>0 



-1/2 /^l 



^at ^ ■ ■ ■ 1 -^at 



a>0 



where the diffusion (x^, pt)t>o solves (★), under the hypotheses ili.). 



3 Asymptotic behaviour of the process (x^) 



a>0 



3.1 An auxiliary function F 

Let us look for a function F = {F^ , . . . , F"^) ^ C'^{W^, W^) such that for I < i < d : 

dF\pt) + dxi - dMi = , (1) 
for some martingale Mf = (M/, . . . , Mf), so that t'^^^ x^t = t"^/^ M*^ - t~^/'^ F^pat). 
Now, Ito's Formula gives : 



dF\pt) 
with 



l<j,k<d 



dt + dMi, 



-1/2 



dMi = [P[1 + v{rt?]) Tin) 9jF'{pt) [dWi + 7?(n) ei dwt]. 

i=i 

Note that, in other words, this means that the so-called infinitesimal generator of the momentum 
diffusion (p^) is 



cj2(r) 



2/5[l + r?(r)2] 



l<j,k<d ^ ^ / j^i I' 



l<j,k<d 

A denoting the usual Euclidian Laplacian of W^. Hence a function F satisfying (1) must solve : 
2B\l7l(r)^] E (^.'^ + ^(^)' ^'■^') 9%F\p) - b{r) j^lP d,F\p) =-p^x f{r). (2) 



Let us take F^ of the form 



F\p) = X ^lJ^{r) = p' X i;pir)/r , 
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and set for r G 



2r6(r) , . , 2r/(r) ^ ^. , 

gir) := ^yV^ H^) ■= l\\ and G{r) 



Then a direct computation shows that 



g{p) dp . 



(3) 



^'p{rt)dWl + 



V/3[l + 77(n)2 

and that Equation (2) is equivalent to : 

d-\ 



n 



■[1 + r/(r)2] 



d-1 



, (4) 



(5) 



Note that, if 6 = / , or equivalently if g = h , then Equation (5) admits the trivial solution 
ipp{r) = r . If d = 1, Equation (5) is easily solved too. But it is not easily solved in the general 
case we are considering, and not even in the case of the diffusion (DH) considered in [DHl] (and 
isotropically extended to higher dimensions) or in [DH2]. 

However, we have the following, whose delicate proof is postponed to Section 5. 

Proposition 1 Under hypotheses iji), Equation (5) admits a solution ipp £ C^(]R+,M) such 
that ipp{0) = 0, and |'0/3('")| = ©(e^'^') , = 0{e^'''') near infinity, for some e' < ^ . 

Moreover, if < f < b , then we have < ip/s < Id . 

3.2 Polar decomposition of the process (p^) and equilibrium distribution u 

Since the diffusion (x^, p^) solves (*), the radial process rt = \pt\ solves : 



with dBt := (1 + r/(n)2)"V2 



Y,0ldWi + r^{rt)dwt 



i=l 



(6) 



As {B,B)t = t, B is a standard real Brownian motion. Consider then the angular process 
= (^^, . . . , ^g) G S'^"^ defined by the time change 9s := Gc-i{s)5 i-e. by pt = x Qct i 



by means of the clock Cj = C{t) :- 



Brownian motion, since it solves : 



/3r| 



ds . The process (©s) G S is a spherical 



del 



1-d 



~Ql\dWl, 



for some standard Brownian motion W = iW^, . . . , W^) G M"'. Hence the infinitesimal generator 
of the diffusion {rt, Qt) = {^t, ©cj is 



a'^ir) 

•= „ ^ o Agd-i , with Cr 

2 pr"^ 



2/3 



d-1 



[1 + i]ir)^] r 



Pg{r) 



dr 



(7) 
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Under this form, it appears that the anisotropy function t] results in a radial drift. 
Set for r G M+ : 



n{r) := exp 



ds 



1 s[l + r]{s)^] 



(8) 



Ll( Si S 

Note that min{r, 1} < ^ufr) < max{r, 1}, and that 0<r<s=^l< —r^ < — . 

/i(r) r 

The radial process (r^) admits the invariant measure v{r)dr , having density on : 

u{r) := a^\r) ^i{rf-^ eT^^^''^ . (9) 

Note that this equilibrium distribution equals the so-called Jiittner one, in the ROUP case 
(we have indeed G(r) = \/l + — 1 in the ROUP and DH cases). 

The hypotheses (H) ensure that v is finite, and then that the radial process (r^) is ergodic. 
Denoting by dQ the uniform probability measure on the sphere S"^"^, and setting : 

7r(dr,(ie) := [j^ v)'^ x v(r)drdQ, (10) 

it is easily seen that vr is an invariant probability measure (or equivalently : equilibrium distribu- 
tion, meaning that the operator A is symmetrical with respect to vr : I ^\ A^2 dix = J ^2 A^i dn 
for any test-functions $1, $2 on M_|_ x S"^"^) for the process (rj , @Ct) — (^ti ©i)? which is then a 
symmetrical ergodic diffusion on M_|_ x S'^~^. 

Lemma 1 For any starting point po = roOo , uniformly with respect to a > , we have : 
r^Epo [|F'(p„t)|2] ^0, ast^oo, forl<i<d. 

Proof Since |F*(p)p = 16** x -0/3 (^')P < V'^l'^) = ©(e^^'^) by Proposition 1, we have 

[|F^(paOH ^ CKroie"^^' , for some constant C. 

Let {Qt) denote the semi-group of the radial diffusion (r^), solution to dtQt = CrQt ■ It is 
known (see for example ([V], chapter 31)) that Qt{ro,r) is a continuous function of {t,r), and 
that (see for example ([V], chapter 32)) Qi{rQ,r) = qi{r)h'{r), for some bounded function qi . 
Hence on one hand we have : 

Er.o[e^^'''1 < sup Qs(e^^'')(ro) < 00 , for < s < 1 , 

0<s<l 

and on the other hand, by the Markov property, for s > 1 we have : 

roo poo 

Er.o[e'^''"1=QiQs-i(e'^'-)(ro)= / ¥.,[e^''-^--]q^{p)v{p)dp < \\qi\\^ e^^'^ u{p) dp < +00 , 

Jo Jo 

by (TC) and since 2e' < (3e . This shows that s 1— > E^ofe^^'''^] is bounded, whence the result, o 

3.3 Asymptotic study of the martingale M 

By Formula (1) and Lemma 1, we are now left with the study of the martingale part (Mt). 
Recall that the coordinates M* of the martingale M are given by Equation (4). 
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3.3.1 Asymptotic independence of the martingales 
Lemma 2 For l<i,l<d,ast^oo we have almost surely : 



lim ^ = 6ii S| , with := 



1 



t^oo t ^' [3d 

Proof The computation of brackets gives easily : 



WpYa'dTT + id-l) / %{'^ + vyid~'(j'd7: 



(5{M\ M')t = 6uSl-{l-6u)Tl' 



with 



SI-.-- 



t ft 

2/„ \ „;./ \2 \ai |2 



aHrs)i^'p{rsf\ei\Us+ / [I + r^{r a\r ^l{r (l 

JO 



3«|2 



) ds, 



a 



[rs) Vl(r-s) [1 + r/(r,)2]-i r~^ - ^['^{r.f 



e\9[ ds 



Setting 
and 



k\r, G) := a'{r) ^'p{ry \e'\' + [1 + v{r)T' a\r) ^^(r) (1 - |0f ) 



f ''(r, 9) := a^{r) [V^§(r) [1 + r]{rfY^ r^^ _ ^'^{j-f^ nial 



and noticing that these functions are vr-integrable by Proposition 1, using Section 3.2 we can 
apply the ergodic theorem, to get the following almost sure convergences : 



lim Si/t 



t— >oo 



k' dn , lim T^'^/t 



ni,l 



dn . 



Now the spherical symmetry with respect to implies that for 1 < i ^ I < d : 

j k'dTT = d~^j |V'^|2(j2ci7r+(l-d-^) j Id-^ a'^d-n =: , and j t'^Tr =0. 

Hence we have got : 



1 /n.ri n,rl\ P-'^: 



t-' {M\ M 



du X S2 = 6u iPd)-^ [vr(aV^l') + (d - 1) 7r(aV|/[(l + 



Consider now the martingale M* defined by : 

M* := (Mi'*, . . . , M^'') := t~'/^ Mat , 
and the Dambis-Dubins-Schwarz Brownian motions i?*'* , such that 

M*'* = B'^\{M''\M''^)a) = B'^\t~^{M\M')at)- 

Applying the asymptotic Knight theorem (see for example ([RY], Theorem (2.3) and Corollary 
(2.4) p. 524-525)), we deduce now from Lemma 2 the asymptotic independence of the martingales 
M* and M\ for 1 < i 7^ / < d, in the following sense. 

Corollary 1 The process (i?^'*, . . . , i?'^'*) converges in law, as t goes to infinity, towards a 
standard d- dimensional Brownian motion B . 
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3.3.2 Convergence of the finite-dimensional marginal laws 



Proposition 2 The finite- dimensional marginal laws of the martingale M* converge, as t goes 
to infinity, to those of the Brownian motion Tip x B, where is the (positive) constant given 
by : 

-1 



"/3 



oo 

dl e-^^^'''^ fi{r)'^-^a{r)-^dr 'x/ il^pir) er ^^^''^ fi^r^-^ h{r) dr . (11) 
Jo 



Recall that tp^ comes from Proposition 1, and that we set in Formulas (3) and (8) : 

r 2r fir) \ r ds 

G{r) = / g{p) dp , h{r) = , p{r) = exp / 

Jo cr^{r) IJi s[l + ri{sy\ 



(12) 



Proof Fix any integer A'^ > 1, positive numbers < ai < . . . < oat , and consider the vector 
random processes : 

X* := ((M''\ M'^^)a, ,BiA , X'=^ := (Tlak, Bi) 

\ 'A. * / l<ij<d,l<fc<7V,s>0 V ^ l<i,j<d,l<k<N ,s>0 

By Section 3.3.1, X* converges in law, as t goes to infinity, to . By the Skorokhod coupling 
theorem (see for example ([K], Theorem (4.30) p. 78)), there exist vector random processes X* 
and X°° satisfying the identities in law : 

(X*) = (X*), 

and such that X* converges almost surely to As a consequence, we get the following conver- 
gence in distribution : 



<d, l<k<N 



B\T}ak)) 



l<i<d.l<k<N 



or equivalently : 



l<i<d,l<k<N P \ ^k) i<i<d,l<k<N 



Note that from Formulas (9), (10), and Lemma 2, we get directly the following expression for 
2 /o°° e-^^W^(r)'^-i dr+{d- 1) il,p{rf e- ^ ^^''^ pirf-^ [1 + r,{rf]-^ r'^ dr 

It remains to derive from this expression the expression (11) of the statement for . This is 
achieved as follows, integrating by parts and using Proposition 1, which implies that 



lim 

1 — >oo 



'ipf3{r)i/j'p{r) e ^^^"^"^ p{r)'^ = , together with Equation (5) : 



ij'Jrf e-^^^'''^ p{rY'^ dr 



dr 



'V'/. (r ) [(^^ - (r )) (r ) - gjji^ 

(5 ijf3{r)e-^^^''^ n{ry-^h{r)dr-{d-l) ijp{r)^ e-^^^'^ p{rf-^ [1 + 7]{rf]-^ r'^ dr . o 
Jo Jo 
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3.3.3 Tightness 

Proposition 3 The family of martingales M* is tight, in C 
of uniform convergence on compact sets of M_|. . 



, endowed with the topology 



Proof Fix any T > 0, and use the Arzela-Ascoli theorem (see for example ([K], Theorem (16.5) 
p. 311)) : the family (r^s M^j, a G [0,r]) is tight, in C ([0,r],R"'), if and only if 



lim limsup E 

t^oo 



-1/2 



sup 

0<a<b<T 
b — a<h 



M„, - M, 



bt 



0. 



Fix i £ {1, . . . , d}, h > 0, and denote by n the integral part of T/h . There exists a standard 
Brownian motion such that (S* being as in Section 3.3) : 



1 



sup |M:,-M,\| = 

fit 0<a<b<T 
b — a<.h 

Setting VF* (u) = sup |VF2|, we have also : 

0<s<u 



1 



0<a<b<T 



sup \W'{Sl,-S:,) 



sup 

ht 0<a<b<T 
b—a<h 



sup - (5^+,, -S:)\< (4,) , 



n<a<T-h 



where 



1 



^ht ■= sup Af^^ , and A^l ■= — 



0<j<n 



SI 



ht ■- V 0+2)hi 



ht 



By the ergodic theorem we have (as in the proof of Lemma 2) the following convergence, as 
t — > oo , valid almost surely and in L^-norm as well : Ajj^ — > 2 S| , which implies the uniform 
integr ability of {Aj^^, ht > 1}. 



Otherwise, by Doob's inequality ^applied to the martingale j l|„<^^j(iM^^ j, we have : 
E 



{Al,)\ < {Al,) ^< 2^E[4J 



whence 



E 



-1/2 



sup \Ml-Ml 



0<a<b<T 
b—a-Ch 

Now, as for fixed h and for any A > 2 S| we have : 



< 2p-'/\ hxE[Al]. 



K[Al] = / P(4, > s)ds < A + / ¥{Al, > s)ds < A + ^ / P(4,^' > s)ds 

Jo JX ._n J\ 



we deduce that 



E[4J < A + 



n 



X Ir.. 



A , as t goes to infinity. 



j=0 



Hence, lim sup 

E[4t] < , and then 



1/9 

lim limsup t ' E 



sup iM^i-M, 



0<a<6<T 
b — a<-h 



^bt\ 



0. o 
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3.3.4 Main result 



Gathering Formula (1), Lemma 1, and Propositions 2 and 3, we get at once the following 
main result of this article. 

Theorem 1 Let (x(,pf) = {xl,pl)i<i<d be a U.'^ x W^-valued diffusion solving the stochastic 
differential system 

{dx\ = f{rt)p\dt 
/ ^ \ -1/2 . . , for l<i<d, 

dpi = - h{rt)p\ dt + a{rt) (/3 [1 + vin?]) W + vin) 01 dwt] 

where W := {W^ , . . . ^W^) denotes a standard d- dimensional Euclidian Brownian motion, w 
denotes a standard real Brownian motion, independent of W, /3 > is an inverse noise or heat 
parameter, and the real functions f, b, a, rj are continuous on and satisfy the following 
hypotheses, for some fixed e > ; 



{TC) a > e on M+ ; g{r) := — ^ ^ > e for large r ; lim e~^ ^ f{r) = for some e' < — 



Then the law of the process (t ^/^X(it)^>Q converges, as t — > oo , to the law of (Y^i^ Ba)^yQ , 
in C(M+,M'^), endowed with the topology of uniform convergence on compact sets of M4. . Here 
B is a standard d-dimensional Brownian motion, and the constant is given by Formula (11). 
This result holds from any starting point (xo,po) (po can also obey the equilibrium law n). 

We deduce now the result conjectured in [DHl], [DH2], and an expression of the limit. 

Corollary 2 Under the same hypotheses as in the above theorem, for any starting point, the 
normalised variance (mean square displacement) t~^ E [|xtp] goes, as t —>■ 00 , towards d x S| . 

Proof By Theorem 1, we have convergence in law of the random variable f~^|xtp, towards 
S^I^Bip. By Formula (1) and Lemma 1, we have only to make sure that for 1 < i < d, the 
following holds : 

t~'K [\Mi\^] = ipty'E[si]^j:}. 

Now, on one hand we already noticed (recall the proof of Lemma 2) that, by ergodicity, we have 
r'Sl ^ . And on .he othe. hand. e™«. .he .n,e . n, .he p.oof 



of Lemma 1 (to show that s E^iJe^^ is bounded), merely using the semi-group (Pt) of the 
diffusion (rt,0cj, solution to dtPt = APt , instead of the radial semi-group (Qt), shows that 
s ^ Ep,j [A;*(p5)] is bounded. Moreover, in the same spirit, by the Markov property and by the 
proof of Lemma 2, we have : 





=-/ 


t 





' Ps{k'){po)ds + I ' Ps{k'){p)ds^qi{p)7r{dp), 



qi being the bounded density of Pi(po, dp) with respect to iT{dp) . It is clear that the first term 
of the right hand side goes to 0. Finally, by the Chacon-Ornstein theorem and by dominated 
convergence, the second term goes indeed to /? S| . o 
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4 Behaviours of , as /? \ and as oo 

Theorem 1 and Corollary 2 show up the interest of the limiting constant d x . Recall then 
from Sections 2 and 3.1 that the processes considered by ([BDRl], [BDR2], [BDR3], [DMR], [DR]) 
and by ([DHl], [DH2]), correspond respectively to : 

{ROUP) h{r) = g{r) = r {1 + r^)-^/^ , G{r) = Vl + ^2 - 1 , r] = , a{rf = 2 , ?/;^(r) = r , 



(DH) h{r) = -— ^ , g{r) = /i(r) = , G{r) = y/l + ^-l , a{rf = 2 y/l+r^, 77(7 

i + r -y/i _)_ 



for some positive (noise or heat) inverse parameter (3 . It is natural to wonder, as in [DHl], how 



behaves the limiting variance T?n , as /3 \ and as j3 00 . 



In the ROUP case, we have simply d x = 2d/ (5 . The variance behaviour is Euclidian. 
In the DH case of [DHl], [DH2], we have by Formula (11) : 



dxTip = 2x 



iPp{r) e-^'^ (1 + r2)-(^+i)/2 r'' dr 



poo 

Jo 

Note that the precise value of ipp is given in Section 5.2 : 

rr poo 

#W=Ci(0/ C2ip)wpipy^h{p)dp + C2ir) Ci{p)wp{pr^h{p)dp, 



with functions Ci;C2;iL>/3 given in Section 5.1. 

In [DHl], for d = 1, after numerical simulations, Dunkel and Hanggi conjecture that 
could be equal to ■ The expression we got above for invalidates this conjecture, and, 
even the asymptotic behaviour near it implies. However, it is true that a non-classical variance 
behaviour occurs. We have indeed the following, whose technical proof is postponed to Section 6. 

Proposition 4 Consider the DH case, for d = 1 , as in [DHl]. Then, we have 

~ 2//3 as (3 Z' 00 ; and, as f3 \0 : ~ ^ — , for some explicit constant A > . 



4.1 Numerical Simulations 



To confirm the validity of our estimates in Proposition 4, invalidating the conjecture of [DHl], 
we performed numerical simulations relating to the DH diffusion, in the case d = 1 . We used the 
Monte-Carlo method, with N = 1000 simulations. For different values of (3 (from 10^^ to 10^), 
we computed x-,(t)-,=i..7v for < t < T = 1000, and then the quantity : 

— 1 ^ 

The following diagram represents our results in logarithmic coordinates. Thus, the horizontal 
axis represents log(l//3), the points * represent the simulated values log(x2 (T)/T) in function of 
log(l//?). 
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The straight hne corresponds to the Euclidian behaviour, the continuous curve to the function 
(3 ^ 2/{j3 + 2), and the dashed curve corresponds to a decrease in log(l//3)~^ for small (3 . 



log(x2(T)/T) 

1 
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These simulations confirm the Euclidian behaviour of the DH diffusion as /3 >> 1 . For small 
/3, the expression conjectured in [DHl] is a good approximation as long as /J > 1/10; however, 
for smaller /? , a divergence appears clearly. On the contrary, the log(l//3)~^-like asymptotic 
behaviour of the limit, which we established above, appears as confirmed. 

4.1.1 The program used for the simulations 

function res=asymp(N,h,D,T) (written in "matlab") 

is the iteration number in the Monte Carlo method. The process x{t) is simulted on [0,T], 
with mesh h. Different values for (5 have been tested. 

Initialisation. Arrays p and x contain the values of p{t) and x{t) for < t <T 

t=0 :h :T; n=length(t) ; r=[] ; 

for k=l :1 :N p=zeros(l,n) ; x=zeros(l,n) ; 

Simulation of Brownian motion 

u=randn(l ,n) ; W=sqrt (2*D*h) *u ; 

Simulation of processes p{t), and x{t) by integration 

for j=l :1 :n-l 

gam=sqrt ( 1+p ( 1 , j ) . *p ( 1 , j ) ) ; 

p(l, j + l)=p(l, j)-(p(l, j))*h + sqrt(gam)*W(l, j) ; 

x(l, j+l)=x(l, j)+(p(l, j)/gam)*h; 

end 

The N simulations of x{t) are placed in the array r 
r= [r ; x] ; end 

Computation of the mean of x^{T), normalised by T 
car=r.*r; limite=mean(car) ; res=limite(n)/T ; 

end of program. 
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5 Proof of Proposition 1 



We are indebted to Reinhard Schafke for this proof, who kindly indicated to us how to proceed 
for Sections 5.1 and 5.2 below. We thank him warmly. Consider first the homogeneous equation 
associated to (5) : 



d-l 



r [1 + r](r)'^] 



(3g{r))c'{r) 



d-l 



C(r) =0. 



(13) 



It has a pole of order 2 at (except for d = 1), and a pole at infinity. Using the fixed point method, 
we construct two solutions Ci and C2 of Equation (13), bounded respectively near infinity and 
near 0. Using these two solutions of the homogeneous equation, a solution ■0^ to (5) is then 
deduced, which vanishes at 0. Finally, we establish the wanted control on tp^,ip'^ . 

I r ds 

Recall from Formula (8) that we set : //(r) = exp 



and min{r, 1} < fj,{r) < max{r, 1}, and 0<r<s=^l< ^^^v < - . 



1 s[l + ri{sy] 



, so that fj, increases 



5.1 Constructing solutions to the homogeneous equation (13) 
5.1.1 Constructing a solution (^i to (13), bounded near 00 

Using hypotheses (TC), fix e > and tq > 1 such that g > e on [ro,oo[. For r > tq , set 



A(r) := 



l-dg/3G{p) 



dp . 



We have 
A(r) < 



-/3e(s-p) 



d-l 



s ds 



dp 



{l + s/p)''-''ds 



p "^dp 



< 



-Pes 



max{l, (1 + s/raf^^} ds = 0{l/r) 



As r ^ 00 , A(r) decreases to 0, so that (up to increase ro) we can suppose that A(ro) < \/{2d). 
On [ro, oo[ , let us define by induction on n G N the functions : (^0 = 1? ^-iid 



(^„+i(r) := 1 + (d - 1) / ^(p)!-'^ e^f?(P) / e-^^(^) pisf-^ ipn{s) [1 + r^{sf]-' s'^ ds 

Jr Up 

We have for r > tq : 

l<'Pn+i{r) < 1 + (d- 1) ||(/9„||lc«[,,(,^oo[ X A(r) < 1 + i ||99„||i«=[^j,_oo[, 
whence 1 < v^n < ||</'n||L°°[ro,oo[ < 2, for any n G N. Then similarly : 

llV'n+l - ^n\\L^[ro,oo[ < WVn - V'n-l || [ro,oo[ X (d - l)A(ro) , 

which allows to apply the fixed point method, to get (i := lim ipn, which satisfies 



dp . 



n — >oo 
L°°\rr.,oo\ 



POO r rco 

l<Ci{r) = l + id-l) fi{p)^-'^e^^^P^ / e-^^^'^ p{s)'^-\iis)[l + 7]{s)^]-h-^ds 

J r Jo 



dp. (14) 
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In particular, as r ^ oo we have : 

Ciir) < 1 + (d - 1) ||Ci||L-[ro,oo[ A(r) 1 , 

hence hm Cii^) = 1 > 

r— >oo 

PCX) 

([{r) = il-d) fiir)^"^ e^^(^) / e~^^^'^ iJ.{sf~^ Ci(s) [1 + visf]^^ s^^ ds < , 

Jr 

then 

Or) + (t,!^ - c;m - = 0. 

This solution can be continued over the whole , yielding d still satisfying (13) and (14) on 

Ml . We have also lim Ci = +oo . 

+ 

5.1.2 Constructing a solution ^2 to (13), bounded near 

For r G [0, 1], set : A(r) := /? / Kpr'^'^ p'^ e'^^^P'^ 

Jo 

We have < A'(r) < /Je^^-^o^'^l / \g{rs)\ sds — > as r ^ , by hypotheses (W). We can 

Jo 

then fix ri G ]0, 1] such that A(ri) < 1/2. On ]0,ri], let us define by induction on n G N the 
functions : (po = 1, and 

(r) :=l + /3 [^{p)""-^ p-^e^'^^P^ 
Jo 

We have G C2(]0, n]), ||(^„+i||ioo]o,r.i] < 1 + A(ri)||(^„||i=o]o,r.j] so that ||(/)„||ioo]o,ri] < 2 , and 

\\4)n+l - (t>n\\L^]0,ri] < {{(pn - (l)n-l\\L'^]0,ri] X Hn) , 

which allows to apply the fixed point method, to get (p := lim cpn , which satisfies : 

L°°]0,ri] 

0(r) = l + P rp^py^-^ g/3G(p) 
Jo 

for any rG]0, ri]. Hence, 

^'(r) = p p{ry^-^ e'^G{r) r-pG{s)^^^-^d+i ^^^^ ^(^^ ^ ^ 0[A'(r)] ^0 as r ^ . 

Jo 

Therefore, (/)(0) = 1, (j)'{0) = 0, and for any r g]0, ri] : 
^"(r) + 

This function can be continued on the whole M+ , into a function cj) satisfying still Equations 
(15) and (16). Set now C2if) ■= r(f){r). It is immediate that ^2 solves (13) on ]R_|_ , and 
satisfies : 

C2(o) = o, C2(o) = i. 



-PG{s) ^^^Y+^^g{s)\sds 



dp . 



-^''^'^p{sf+^g{s)<Pn{s)sds 



dp . 



-^^^''^ p{sf+^ g{s)(l){s)sds 



dp=l + 0[K{r)] (15) 



-r— — —-^ + - - (3g{r) 
r [1 + ri[ry\ r 



cP'{r) - P^<P{r) =0. 



(16) 
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5.1.3 The Wronskian Wj3 of Ci)C2 
Consider the Wronskian : 

Wf3 := C1C2 - C1C2 on 

We have 



W'f, = C1C2 - CfC2 ={P9- , r, ) X "'/^ 



d-l 
[1 + 772] /d^ 

so that 

W/3(r) = a/s fi{r)^~'^ e^^^^""^ for any r > 
and for some constant . As Ci ^ 1 > C2 > > Ci < > C2 > near , we must have a/3 > . 



5.2 Constructing a solution tpis to Equation (5) on 



For any continuous function k on M4. , such that Hm e ^ ^k{r) = for some e' < ^ , and 
for any < r < 00 , set : 



^{k){r):=Ci{r) Q2{p)wp{p)-^ k{p) dp + C2{r) Ci{p) wp{p)-^ k{p) dp . 

Jo Jr 

Note that ^'(A;) is well defined; we have indeed, using that g > e on [ro,oo[ : 

Up)wfs{p)-'\k{p)\dp = 0{l) / p{p)'-''e-^'"^PK^'P/'dp <oo. 



' ro -J ro 

Note also that by (TC), we can take in particular k = f3h (recall Formula (3) defining h). 
Moreover, again for < r < co , we have : 



00 



^iknr)=C[{r) I C2{p)w^{pr'k{p)dp + C',{r) I Ci{p)wp{p)-^k{p) dp , 

and 

^(^)"W=CrW / Mp)wp{p)~^k{p)dp^C2{r) I Ci{p)MPr'Hp)dp-k{r), 







r POO 







SO that ^(k) solves on Equation (5), with k instead oi [ih. Near 0, we have C2('') ~ i" ■ 
Otherwise, noticing that (C1/C2)' = ~W[3/C,2 > have for r > : 

/' ^ = ^(0 - ^(1), Ci(r) = ^(1) Ur) + Ur) f ^ ■ 

Jr C2 "^2 (,2 C,2 Jr C2 

Hence, near we have : Ci('") ~ ^^.p / p{s)^~'^s^'^ds < ap p{r)^^'^, and then "^{k){r) = 0{r). 

Jr 

In particular, we have ^'(A;)(0) = 0. Using Cii^) = (^(^pi^y ^^^ the expression of (^[ (recall 
Section 5.1.1), we get at once ICK'^)! = Oip{r)^'''^/r] near 0. We have therefore near : 



\C'i{r)\ r C2{p)wp{pr'k{p)dp = 0(^i(r)i-7r) f\{pf-^ pdp = 0{r) 
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and 

C2(0 / Ci{p)wp{p)-^k{p)dp = o{i), 



POO 

whence ^(A:)'(0) = / Ci{p) Wfsipy'^Kp) dp £ R . Using (5), we get also <(0) G M . 
Jo 

Setting ipij := (3'${h), we have thus V/3 S C^(M+), tppiO) = 0, and ipp solves (5) on 

5.3 Estimates for tp/^ and tp'^ near oc 
Recall from Section 5.2 that 

pT poo 

^p{r)=(3Ci{r) C2{p)wf3ipr'h{p)dp + (3C2{r) Qi{p) wp{p)-^h{p) dp 



and 

4(r)=/?C((r) / Up)^^p{pr^Hp)dp + [iQ'2{r) I Ci{p)w^{p)-^h{p) dp . 







Near infinity, we have on one hand : Ci ~ 1 ^'^'^ (by Section 5.1.1) 

\Q'i{r)\ = o{p{rf^'^'^j e-'^''^'-'''^ p.{sf-^ s-'^ ds = Oir^"^) j e-^^' [1 + s/r)'^~Us = 0{r^^). 

On the other hand, by (Ji.) we have also near infinity : 

2 , , , / 2 , , 2(d-l) \ 

Mr) - ^ w,{r) =[l-- gir) + j Mr) < , 

r 2 

so that wp < — iwf3{r) - Wfj^ro)) = C'(u;^(r)). 

Jro P ^ 

Noticing that (C2/C1)' = wp/Cf , we have for any r > : C2{r) = ^(1) Ci{r) + Ci{r) [ ^ . 

u Ji Cl 

This implies that C2(^) = near infinity. Therefore, by (TC) there exists an e' < e/3/2 

such that : 

r C2(p) «^Mp)-'^(p) dp = f 0{e^'P) dp = 0{e^'^) . 
Jo Jo 

To control the second integral in the expressions for ipis and ip'^ , observe similarly that : 

poo poo 

wp{r) / Ci(p) wp{p)-'h{p) dp = p{rf-^ / e-^[^(^)-^«] p{p)''-' 0{e''P) dp 



Oie''"^) p{r)^-^ / e-'^^P-"^^/^ p{p)'^~^ dp = 0{ef''') / er'^P/^ (1 + p/r)'^-^ dp = Oie^'"") 





Then by definition of , we have (^2 = + ~~ Ci 1 whence C2 ~ '^P near infinity. 

Ci Ci 

As a conclusion, gathering the above we have indeed, for some e' < e/3/2 : |^/;^(r)| = ©(e"^ 
and I "0^(^)1 = 0{e^ ''') ^ for large r . 
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5.4 We have p ^{g) = Id on R+ , and 0<^f3<Id if < / < 6 

Set "0/3 •= P^is)^ s-iid note that the identical function Id solves (5) if g = h . Hence, the 
function r tjjp{r) — r solves the homogeneous equation (13), so that 

tpfsir) — r = c C2(^) + c' Ci{r) , for some real constants c, c and for any r > . 

/r 

must have c > . Otherwise, integrating by parts, near infinity we have : 

J r J r 

SO that 

roo poo 

(3 / Cl{p)^^{pf-^e~P''^P^g{p)dp ^(3 / p{pf-^ e'^^'^P^ g{p) dp ^ pirf-^ e-^'^'^''^ 

J r J r 

Hence, by definition of ipi^ and ^ and by the above, we have near infinity : 



^l^p{r) = 0{l) g{p)dp + 0{wp{r))ii{r 



whence 



ce 



PG{r)/2 



o[cC2{r)) = o['4)p{r) -r) = o[e 



Therefore V'/sC'^) = r on M+ , as wanted. Finally, if < / < 6 , then Q < h < g ,so that ^'(/i) > 
and ^{g-h) > by Section 5.2 above, and then by linearity of ^ : Id = (3^{g) > f3^{h) = tpp . 
o 



6 Proof of Proposition 4 

pr r pc 

For d = 1 (and any 77), it is immediate from (5) that 'il)p{r) = [3 I I 

Jo Up 

so that in the DH case, the limit expresses as : 



dp, 



S2 =2/? 



„/3(l-vTf^^) 



■ (ix X 



l + y2 



■dy 



where 



dz, Kp := 2(3 I Ip{xfe^^^^-^)dx, Ip{x) := 



ye 



■dy . 
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6.1 Behaviour as f3 ^ oo 

Integrating by parts yields : 



1 yy<e 



-/3(\/T+y2-i) 



13 J.. (l + y2)3/2 



■dy . 



As 



dy 



J X 



1 + X2 I J,^ (l + y2)3/2 



(l + y2)3/2 



we get 



Hence, 



H.fe^i^-^) = '-^r^^^ X [1 + O (1//?)] . 



2 /■ p[i-VT+^) 

— dx X [1 + O (1//3)] . 



Setting u = P [ VT+~x^ — 1 j , we get : 

2 ^ g/3(l-v^T4^) 



e ^ du 



PJr+ 1 + (n + /J)v^u(n + 2/3) 

by dominated convergence, as (3 ^ oo . We have similarly : 



^ —3/2 



/3" 



/3 



„/3(l-VT+F) 



e du 



yju{u + 2/3) 



^ --1/2 



whence 



e2 = J-^i^^ ~ 2/(3 . 



6.2 Behaviour as /? — 

We have 



J/ 



/3 



1 e-^"'dn /-i/v^ e-'/''dt 



V^u('u + 2/3) Jr+ Vii2 + 2 
where we performed the change of variable /Ju^ = l/t^ . Hence, as /? — > , we have : 

"i/v^ (e-i/t^ _ 1) dt 



J« = 2 



(it 



1 t Vl + 2/?t2 
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+ Ci + o(l), 



for a positive constant Ci . Integrating by parts yields then : 

l/v^ 

Setting u = t ■, we get : 



+ 2/3 



J 1 



tlogtdt 



1 (1 + 2/3*2) 



3/2 



log 



1 tyrrw kvp 



1 



+ 2 



1 



:l0g 



1 



V3 "VV^ 

Now, as /3 ^ we have 
u du 



+ 2 log 



du 



^ u\og{u / y/]3) du 

VP (1 + 2^x2)3/2 

^ ti log ti dii 



V v^y y v/j (1 + 2^/2)3/2 (1 + 2^2) 



+ 2 



V/I (1 + 2n2)3/2 



and 



Hence, 



1 

2Vl + 2'u2 
^ u log n dn 

v//7 (1 + 2tz2)3/2 " 
l/V/? fit 



^" 2 V7r+1^"7! 



\3/2 ■ 



V3 



Ca :=2 



^ u log u du 



h t Vl + 2/3*2 
and J/3 ~ 2 1og(l//?). Otherwise, 



/o (l + 2n2)=^/2 
log(l/y^) X [l + o(l)]+C2, 



K^ = 2p[ Ip{x)'e^(^-')dx = 2P [ (r e^^'-^^-Ye^i-^-')da 

Jr+ Jr+ \Jvt+^ " / 



oo / roo 



2/3 

'1 \Jy 

Observing that 



dy = 2/3 



oo / />oo 



1 \Jy 



3u ^\ ^ y i^^/ 





2 


ry/f^ ^ du 


2 


= 2 








+ 


' r,-pudu 


Jy U 




4 n_ 






Jy n_ 











1 \ e-?' e"2y 



and then that 



log 



-I3u 



du 



y dy 



log 



/3 



-1 



y/f^ a du 



y 



y dy 

77^ 



+ 0(1) 



changing /3y into y we get : 



(l„gi)"K, = 2/(logi)"A, + 0(/3), with A„:=I^(^J^ 



°° /"^/^ -ud^\ ye^dy 
u 
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Finally, setting x = y^y^ _ ^2 ^ -^^g have by dominated convergence : 

Jo \Jyp^W ^ / Jo \Jx U J 

whence 

e2 = J-^Kp ~ yiog(l//3). o 
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